Draft version June 10, 2010 

Preprint typeset using 1^1^^ style emulateapj v. 11/26/04 



CRITERIA FOR FLUX ROPE ERUPTION: NON EQUILIBRIUM VERSUS TORUS INSTABILITY 



o 

(N 

a 
6 



p. Demoulin and G. Aulanier 

LESIA, Obscrvatoirc dc Paris, CNRS, 5 place Jules Janssen, 92190 Mcudon, France 
Draft version June 10, 2010 

ABSTRACT 

The coronal magnetic configuration of an active region typically evolves quietly during few days 
before becoming suddenly eruptive and launching a coronal mass ejection (CME). The precise origin 
of the eruption is still debated. Among several mechanisms, it has been proposed that a loss of 
equilibrium, or an ideal magneto-hydrodynamic (MHD) instability such as the torus instability, could 
be responsible for the sudden eruptivity. Distinct approaches have also been formulated for limit 
cases having circular or translation symmetry. We revisit the previous theoretical approaches, setting 
them in the same analytical framework. The coronal field results from the contribution of a non- 
neutralized current channel added to a background magnetic field, which in our model is the potential 
field generated by two photospheric flux concentrations. The evolution on short Alfvenic time scale 
is governed by ideal MHD. We show analytically first that the loss of equilibrium and the stability 
analysis are two different views of the same physical mechanism. Second, we identify that the same 
physics is involved in the instability of circular and straight current channels. Indeed, they are just 
two particular limiting case of more general current paths. A global instability of the magnetic 
configuration is present when the current channel is located at a coronal height, h, large enough so 
that the decay index of the potential field, 9 In |Bp|/91nft, is larger than a critical value. At the limit 
of very thin current channels, previous analysis found a critical decay index of 1.5 and 1 for circular 
and straight current channels, respectively. However, with current channels being deformable and as 
thick as expected in the corona, we show that this critical index has similar values for circular and 
straight current channels, typically in the range [1.1,1.3]. 

Subject headings: Sun: corona — Sun: filaments — Sun: flares — Sun:magnetic fields — Sun: photo- 
sphere 
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1. INTRODUCTION 

A coronal mass ejections (CME) is the consequence of 
the sudden dcstabilization of a part of the coronal mag- 
netic field. The eruption is preceded by a long phase 
(days to week) during which the magnetic field is pro- 
gressively stressed and free magnetic energy builds up. 
The configuration typically grows quasi-statically (with 
velocities well below the Alfvcn velocity). At a point of 
the evolution, in a few minutes up to an hour, the sys- 
tem becomes very dynamic, with a global upward mo- 
tion, as traced by the evolution of the cold plasma in 
the associated filament and of the hot plasma in coronal 
loops. Later on, a significant release of magnetic energy 
occurs, and a flare is typically observed. If the down- 
ward magnetic tension of the covering magnetic arcade 
is weak enough, the erupting plasma and magnetic fleld 
is launched towards the interplanetary space as a CME. 
In summary, the CME phenomena occurs in four main 
phases: build-up, instability, accele ration, and propaga - 
tion They have been reviewed by iForbes et all ()2006f ) 
and lVrsnakI (l2008l ). 

The two last phases are the most spectacular ones, 
so they are better constrained by observations and they 
are more deep ly modeled, in part i cular with MHD sim- 
ulations (e.g! lAmari et al.l 120031 : iFan fc G\hsa^ 120071 : 
iTorok fc Kliemll2007f) . The first phase is a slow evolu- 
tion and it is usually difficult to characterize in obser- 
vations what arc the generic key points which lead to 
eruption. The main physics which emerges from obser- 
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vations is the presence of new emerging magnetic flux, 
progressive dispersion of the whole flux, the build up of 
a very sheared field in the vicinity of the photospheric 
inversion l ine (PIL), and the cancelatio n of flux at the 
PIL (e.g. &an Driel-Gesztelvil [i99& Green et al.l 120021 : 
Ivan Driel-Gesztelvi et al.ll2003f ). At the least the three 
last are physically related since flux dispersion lead to the 
convergent flows towards the PIL, increasing the mag- 
netic shear and forcing flux cancelation. This also implies 
the build-up of a flux rope with J-shaped coronal loops 
transformed by re c onnection into S-sha.ped loops (e.g. 
Moore et al.l [l995l : iGibson et all l20M iGreen fc Klieinl 



20091 ). 



A still open issue is why does the magnetic configu- 
ration erupt? There is usually no evidence of a large 
amount of new magnetic flux (with a magnitude com- 
parable to the pre-eruptive flux), so the eruption is 
not driven by the sub-photospheric evolution but rather 
the coronal magnetic conflguration becomes unstable 
at some point during the slow evolution. During this 
phase, magnetic reconnection is probably involved as a 
key mechanism for the progressive transformation of the 
magnetic conflguration. However, an ideal instability is 
thought to initiate the CME since the upward accel- 
eration phase starts before the imp ulsive phase of the 
flare in the majorit y of events (e.g. I K abler et al.l 119881 : 
iMaricic et 311120071 ). Later on, magnetic reconnection 
plays a key role in the eruption as the peak of the up- 
ward acceleration is typically found correlated with the 
peak of the hard X-rays and of the time derivative of soft 
X-rays flux (e.g. iNeupert et al]|2001l : IZhang et all 120011: 
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Fig. 1. — (a) Schema of a current channel in the corona and its 
image current below the photosphere (located at 2 = 0). The cur- 
rent channel has a radius a. (b,c) Particular cases with a circular 
and straight current channels, respectively. 



IVrsnak et ani2004t iTemmer et al.ll2008| ). 

Magnetic reconncction also plays a key role during the 
first phase as it permits the progressive transformation 
of very sheared field lines int o a twisted flux rope. How- 
ever, the MHD simulations of lAulanier et al.l ()2010l) have 
shown that magnetic reconnection at the photospheric 
level, or later on in the corona below the flux rope, is not 
directly responsible for the onset of the eruption. The 
conflguration rather gets unstable when the flux rope 
reaches a height where the potential field, associated to 
the distribution of the photospheric magnetic flux, is de- 
creasing fast enough with height. This relates the onset 
of the eruption in MHD simulations to a series of analyt- 
ical studies, as summarized below. 

The equilibrium of a flux rope within a c orona l field 
was first considered by iKuperus fc Raadul ()1974| ) with 
the following simplifications. In Cartesian coordinates, 
they modeled a flux rope with the magnetic field cre- 
ated by a straight line current of intensity /, located at 
a height, z = h, above the photosphere located at z = 
(Figure [Tfc)). The flux rope field is added to a simple 
background magnetic field: a potential field, Bp, asso- 
ciated to a bipolar photospheric magnetic field. They 
included the observed insignificant evolution of the ver- 
tical component, B^, at the photospheric level during the 
instability phase of a CME by introducing an image cur- 
rent of intensity — /, located at the height z = —h. The 
physical result is that two oppositely directed Laplace 
forces are acting on the coronal line current: one from 
the potential field. Bp, and the other from the magnetic 
field created from the image current (or the equivalent 
surface current at z = 0). The equilibrium is then de- 
s cribed by a curve Ijh) . 

Ivan Tend fc KuperusI ()1978f ) showed that, I{h) is an 
increasing function of h at low height and that it has 
a local maximum if the horizontal component of Bp or- 
thogonal to the line current, |i?p,x|, decreases fast enough 
with height. Then, supposing that the current / can 
be increased progressively to larger values during the 
buildup phase, a loss of equilibrium occurs at ft, = ftcrit. 
defined by the maximum of the function I{h). This 



occurs where |i3p,x| decreases faster than 1/h. This 
model was later developed within a circuit theory, in- 
troducing an electric potential and a resistance in the 
circuit to describe the temporal evolution of /. The 
current sheet which forms below the erupting flux rope 
was described by another electric circuit. This pro- 
vides a set of coupled equatio ns which describes the 
main phases of an eruptive flare (iMartens fc Kuinlll989l : 
Ivan Ballegooiien fc Martens! Il989f )^ Next, the^ back- 
ground potential field, Bp, was replaced by a linear 
force- free field to in clude, as observed, a sheared field 
(jAmari fc Alvl Il989f ) . A loss of of equilibrium is also 
present if |i?p.y| is decre asing fast enough with height 
i|Demouhn fc Priest) [l988l ). 

In the above theory, it is supposed that the current in- 
tensity, /, is the main driver of the evolution and that it 
can be increased to arbitrary large values. However, the 
coronal physics is not precisely described by any circuit 
theory because the electric current is a consequence of 
the force bala nce together with MHD constraints (e.g. 
lParker|[T996ah . The effect of magnetic reconnection on 
the system (through changing field line connectivity and 
energy release) before the flare/CME acceleration phase 
is minimal and therefore the stability of an equilibrium 
is typically tested in ideal MHD. The conservation of 
the coronal magnetic flux passing below the flux rope is 
typic ally used to set a constra int on the current evolu- 
tion. lAnzer fc Ballesterl (|199Cl[ ) claimed that there is no 
longer a los s of equilibrium with this ideal MHD con- 
straint, but iDemoulin et al.l ()1991l ) found that this con- 
straint mostly displaced the loss of equilibrium point to 
a larger height (slightly after the maximum of t he 1(h) 
curve) , in agreement with the MHD simulation of lForbea 
(|2000( l. The following developments have shown that a 
loss of equilibrium is typically present if the flux rope 
radius is thin enough. It occurs, for example, when the 
photospheric polarities are subject to converging motions 
toward the PIL, or when their magnetic flux is decreased, 
even when ideal MHD is assumed during the full build- 
up phase, ll senbcr g et al.lfl993t iForbes fc Isenber^ll991l : 
iForijes fc Pr iest 1991). 

When the current channel is curve d, an extra forc e 
is present, called the hoop force (e.g. iBatemaril Il978| ). 
The electric current of a curved channel creates a mag- 
netic fleld component orthogonal to the channel. This 
implies an outward Laplace force (away from the curva- 
ture center). In terms of the magnetic field, this force 
is due to the over magnetic pressure of the azimuthal 
field component present in the direction of the curvature 
center. This for ce is a t the heart of the magnetically 
driven mo del of IChenI (Il989f) and subsequent de velop- 
ments (e.g. lGarren fc Chenlll994l : lKraU et al.][2000l) . This 
force is also present in MHD models where the straight 
line curr ent of previous paragraph is replaced by a ring o f 
current (iLin et al.lll9Ml2002l : lTitov fc Demoulinlll999t) . 
As for the above cartesian models, an ideal MHD evo- 
lution during the build-up phase also typically, but not 
always, lead to a loss of equilibrium. The differences 
will be discussed in Sections |3| and |31 A non equilibrium 
point could also be present when line-tied conditions are 
imposed at the photosph eric footpoints of the flux rope 
(jlsenberg fc Forl3esll2007t l. 

From another point of view, iKhem fc Tor"6l3 (|2006h 
studied the stability of a toroidal current ring immerged 
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in a background po tential fie l d, Bp . Extending the re- 
sults summarized in lBatemanl ()1978[ ). they derive that an 
instability occur when the background field component 
orthogonal to the torus (|-Bp_x|) decreases faster than 
1 / h^/"^ with a correction dependi ng on the torus aspect 
ratio (major over minor radius). iKliem fc Torokl ()2006l ) 
called it the "torus instability". They analyzed cases 
where the electric current / was held constant or fixed 
by the conservation of the total magnetic flux within the 
torus hole. 

The aims of present paper is to revisit the above studies 
to analyze their relationships. Is this "torus instability" 
different from the loss of equilibrium found in previous 
studies with a toroidal current? H as it a different phys- 
ical o rigin than the instability of Ivan Tend fc KuperusI 
(|1978f ) obtained with a straight current channel? The 
Cartesian and axisymmetric models are first written with 
the same formalism in Section We compare the loss 
of equilibrium and the instability approaches in both ge- 
ometries in Section [3l before comparing the criteria of 
instabilities taking into account a finite current-channel 
width which evolves during the perturbation ( Section |4]). 
Finally, in Section [5l we summarize our results and dis- 
cuss the relationship between the Cartesian and axisym- 
metric models, as well as between the loss of equilibrium 
and the torus instability. 

2. BASIC CONCEPTS 

The coronal magnetic field B can always be is written 
as the sum of a background magnetic field Bp and the 
magnetic field Bj, created by localized coronal currents 
and their images below the photosphere (Section 12. 2|) . 
Bi has a vanishing component at z = by con- 
struction. All the analytical models use this field de- 
composition with various approximations for Bj. Bp is 
most frequently taken as the unique potential field as- 
sociated with the photospheric "normal magnetogram" 
Bz(x,y,0) and with vanishing strength at infinite dis- 
tance. We consider this case below, however, it is worth 
noting that, in general. Bp can also incorporate dis- 
tributed coronal currents (e.g. Bp could be a linear forcc- 
free field, and Bj needs t o be self consistently computed, 
iDemoulin fc Priest! 1 19881 ). 

2.1. Modelization with Concentrated Currents 

In order to have a set of equations solvable analyti- 
cally, simplifications in the magnetic configuration need 
to be made. A first one is to suppose that the elec- 
tric current is restricted to a thin channel in the corona 
(Figure [T] (a)). More precisely that its typical thickness, 
2a, is small compared to the spatial scales of Bp and 
to the local radius of curvature of the current channel 
axis. In an active region, the magnetic shear is typically 
concentrated around the PIL, while the surrounding ar- 
cade is more potential, so electric current are stronger 
around and above the PIL. With reconnection of sheared 
loops at the PIL, an importa nt fraction of the cur rents 
are inside a twisted flux tube (jAulanier et al.l[201ClL and 
references therein). Then, the introduction of a concen- 
trated current channel is motivated by observations and 
MHD simulations, however it is still an important sim- 
plification (e.g., neglecting the effect of more distributed 
currents as well a the presence of narrow current layers 
and sheets). 



The approximation of a thin current channel al- 
lows to separate approximately the magn etic equilib- 
rium in a internal an d external equilibrium (|Chenlll989l : 
llsenberg et al.|[T993[) . This splitting is better achieved 
as the current channel is thinner. For the external equi- 
librium, the Laplace force, integrated over the channel 
cross section, vanishes, so that there is no average mag- 
netic field component orthogonal to the current channel. 
The internal equilibrium is solved locally, in the chan- 
nel cross section. A twisted flux tube has generically 
both toroidal (axial) and poloidal (azimuthal) magnetic 
field and electric current components. The equilibrium 
is typically solved in cylindrical coordinates with a bal- 
ance between the total magnetic pressure gradient and 
the tension of the poloidal magnetic field (force-free field 
solution). This internal equilibrium is not t he object 
of pre sent paper, and we refer to the work of iLin et al.l 
(fT998l) . 

Another important simplification is the absence of a 
neutralization, or return, current around the direct cur- 
rent. Such return current, of opposite direction and 
with the same magnitude as the direct current, is ex- 
pected to be present in a magnetic field formed by emer- 
gence or induced by loc alized boundary (photospheric) 
motions ()Parkeij[l996bl ). In both cases the complete 
neutralization is due to a vanishing circulation of B 
around a large path enclosing the current channel. Re- 
turn currents are indeed present in MHD simulations, but 
when a significant magnetic field component is present 
along the PIL, the direct curre nt has a lar ger mag- 
nitude than the ret urn current (jTorok fc Klic ni ,20()J; 
lAulanier et al.|[2005l) . Indeed, only a partial neutraliza- 
tion was typically reported in sun spots ()Wheatlandll200Cll : 
IVenkatakrishnan fc Tiwarill2009l ). With partial neutral- 
ization, the current intensity, /, is the non- vanishing sum 
of the two opposite currents. 

We argue that the occurence of non fully neutral- 
ized currents should be a common feature in solar ac- 
tive regions. The contrary would imply that current- 
carrying flux tubes should be fully surrounded by po- 
tential fields, not only high up in the corona, but also 
low down around PILs. It is worth noticing, however, 
that some MHD models for solar eruptions clearly do 
not require a net current (the magnetic breakout, tether 
cutting and flux disappearance m odels, respectively ad- 
dressed bv lAntio chos ct al. 199^; iMoore fc RoumeliotisI 
Il992t lAmarreral.ll2000[ ). as they are based on a grad- 
ual diminishing of the tension of the background field, 
irrespectively of the distribution of electric currents at 
lower a ltitudes. Many other models exist (as r eviewed 
e.g. in lForbes et"al] 120061: lAulanier et al]|2010[) . but it 
is conceivable that several actually fall into the physical 
frame studied in the present paper. 

2.2. Image Current 

The lack of significant photospheric magnetic flux evo- 
lution during the initiation of a CME can be modeled 
with the introduction of image curren ts below the photo- 
sphere . The straight channel case of iKuperus fc Raadul 
(fl97l can be generalized to any channel shape as il- 
lustrated in Figure HJa). Let the current vector be 
{Ix, ly^Iz) at a generic point (x, y, z) of the corona, then 
the introduction of the image current {—Ix,—Iy,Iz) at 
the image position [x, y, — z) implies that the vertical 
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component of the magnetic field at the photosphere, 
Bz{x,y,0), is unaffected by the presence of the current 
channel. 

The introduction of this image current is a particu- 
lar case of the technique of images to impose a particu- 
lar bound ary condition i n the setting of electrodynamic 
problems (|Jacksonlll975l chap. 2 and 5). Physically, the 
coronal current path is closing in a complex set of hori- 
zontal photospheric currents, which create the same coro- 
nal magnetic field than the image current. 

2.3. Magnetic flux 

We suppose that photospheric evolution and magnetic 
reconnection are negligible during the instability phase, 
then the total coronal magnetic flux passing below the 
current channel is conserved. This flux is the sum of the 
flux of Bp and of Bj. By the symmetric construction of 
the image current, the coronal flux of Bi is half the flux 
enclosed by the full current channel (coronal and image 
current), so it is equal to LI/2, where L is the externa l 
inductance of the full current channel (e.g. lJacksonllQTSl ). 

The total inductance of a circular channel of main ra- 
dius h and of small radius a (Figure [TJb)) is 



(1) 



where /io is the magnetic permeability, and l\ is the nor- 
malized inter nal inductance , per unit length, of the cur- 
rent channel (jGroveilll946[ ). h take the value of 0, 0.5 
or 1 for a current concentrated at the border of the 
torus, uniformly distributed within the cros s section, 
OTwith a linear force-free field equilibrium (|Lin et al.l 
Il998l and references therein). Equation ([T]) is simple, 
but still a good approximation of more complete ex- 
pressions. For e xample, it is close to the series expan- 
sion tabulated by iMalmberg fc RosenbluthI (|1965l) in the 
range 0.1 < a/h < 1 for li = 0. M oreover, the expres- 
sion with elliptic integrals given by iRamo et al.l (|1994 
p. 193) is also close to Equation ([T]) with /i = in the 
range 0.15 < a/h < 1, while for lower a/h values it is 
closer to the case = 0.5. Finally, for the coronal mag- 
netic flux passing below the current channel, the external 
inductance is required, so the flux is LcI/2 with li = 0. 

With the same notations (Figure [ijc)), the total induc- 
tance of a straight channel and its image, for a length Ay 
along the channel, is 



(2) 



where l^ has the same value than in the above circular 
case. The main difference with the circular case, is that 
Ls is only weakly dependent of h so that Lg is almost con- 
stant during a global perturbation of the channel (modi- 
fying the height h) . As above, the coronal magnetic flux 
passing below the current channel is L^I /2 with li = 0. 

2.4. Magnetic Self Force 

Prom the Biot and Savart law, a current channel gen- 
erates a magnetic field at any point of the space, in 
particular in the current channel. The component of 
this field orthogonal to the current channel induces a 
Laplace force. The direct calculation of this force from 
the Biot and Savart law is very cumbersome, even for 



a simple torus geometry. In practice, this force is com- 
puted by equaling the work of this force to the change of 
the magnetic energy (LJ^/2) during an elementary dis- 
placement, preserving the magnetic flux encircled by the 
current channel, LI, so that there is no inductive effect 
(iShafranovl 119661 : [G arren fc Chenl[l99l . For a circular 
current channel, this implies an outward radial force fc, 
per unit length along the channel, given by: 



fc 



PdLc 

2 dh 



j2 A^o 



In 
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(3) 



The force fc is the Laplace force between the toroidal 
(axial) current and the poloidal magnetic fleld. There is 
also the Laplace force between the poloidal current and 
the toroidal magnetic field inside the twisted flux tube. 
Taking into account the internal equilibrium, in the limit 
a/h « 1, this second force only gives a small correction 
to the previous force, as the total force is obtained from 
Equation Q by replacing li by l[ = k — Sh with Sh ~ 1 
(|ShafranovE966l) . Then, we write the total force as rcl 



with 



inh 



(4) 



The local outward force rcl^, called hoop force, has its 
origin in the magnetic fleld created by each of the ele- 
mentary part of the circular ring, with the largest con- 
tribution coming from the closest currents. Indeed, it 
has a logarithm divergence as the small radius, a, be- 
comes smaller. For a « h the logarithm term slightly 
dominates in Equation (l4l). If the cur rent channel is 
fully in the corona (e.g. iLin et al.l [19981) . Equation Q 
includes only the self-force of the coronal current chan- 
nel, whereas if the current channel is only half in the 
corona fe.g. iLin et al.ll2002l Figure mb)) Equation (gj) 
also includes th e force from the image current. Finally, 
iGarren fc ChenI (|1994D have generalized these results to 
current channels with arbitrary shapes. 

For a straight current channel, located at the height h 
above the photosphere, the magnetic force can be com- 
puted as above (Equation ([3])) with Lc replaced by Ls or 
the magnetic field of the image current can be computed 
first by the Biot and Savart law. The repulsion function 



Mo 
47r/i 



(5) 



has indeed a similar form as Tc (Equation (j4|)) since 
the logarithm term provides only a weak dependance on 
h/a. 

3. LOSS OF EQUILIBRIUM AND INSTABILITY 

3.1. Magnetic Field Evolution 

The MHD evolution of the coronal magnetic fleld of 
an active region, outside flare times, has typically three 
time-scales, as follows. The shortest time-scale, ta, 
called the Alfvcn time, is given by the typical time that 
Alfven waves require to cross the coronal field configura- 
tion. Outside of the eruption period, the magnetic con- 
figuration is not significantly changing on such time scale 
and Alfvcn waves transport the magnetic stresses (shear, 
twist). On an intermediate time-scale, tb, the coro- 
nal magnetic stress is significantly changing, for example 
because of sub-photosphcric torsional Alfven waves are 
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bringing twist to the coronal field, or because of recon- 
nection between sheared magnetic loops is forced by con- 
verging flows at the PIL. On this intermediate time-scale, 
coronal currents significantly change, but not the global 
photospheric distribution of the magnetic field compo- 
nent normal to the photosphere (i.e. the "normal magnc- 
togram"). On the longest time-scale, re, both the coro- 
nal currents and the "normal magnetogram" are evolv- 
ing. 

For a mature AR, with a magnetic fiux of w 10^^ Mx, 
the above time-scales are typically of the order of ta ~ 
few minutes, tb « few days, and a tc « week to months. 
They are well separated, a useful property for an ana- 
lytical study since it allows to isolate the main physics 
involved for each time-scale. Thus, we can use different 
approximations to analyze the magnetic field evolution 
on these time-scales, as follows. The global stability of 
the equilibrium can be tested with ideal MHD on the 
shortest time-scale ta- On the intermediate time-scale 
Tb , the build up of coronal current can be studied within 
a given potential Bp (computed uniquely from the "nor- 
mal magnetogram" , using appropriate boundary condi- 
tions on the sides of the domain). Finally the evolution 
on the longest time-scale tc is mostly studied with MHD 
simulations or with analytical studies imposing an ideal- 
MHD evolution. 

3.2. Loss of Equilibrium with a Circular Current 
Channel 

First, we re vi sit the model proposed by 
iTitov fc DemoulinI ()1999( ). keeping it the simplest 
as possible. The potential field Bp is created by two 
magnetic sources, of equal flux but of opposite sign, 
located at a; = ±i?, y ^ z = and by a line current 
located at ?/ = z = 0. The magnetic field of this line 
current adds a contribution to the coronal potential 
field, but it does not modify the equilibrium, nor the 
equilibrium perturbation studied below, so this field 
is not present in the following equations (the aim of 
this magnetic field is only to have a finite twist in the 
coronal field). A torus of electric current, centered 
atx = y = 2: = 0, is introduced in the plane y-z 
(Figure [TJb)). Since the photospheric boundary is set at 
z = 0, the half part of the torus at z < represents the 
image current. 

The magnetic configuration considered is axisymmet- 
ric (around the x-axis), then the magnetic force on the 
current channel is radial in the y-z plane with the same 
value along the current channel: 

/ = re/2 + Bp.,/, (6) 

where rc is given by Equation The x-component of 
Bp at the current location is; 

Bp,^ = -4(PD{h^ + D^)-^^\ (7) 

The equilibrium current, /cq. is given by / = in Equa- 
tion ®: 



/cq.(/l) = 



167r (l)Dh{h^ + D^)-^/^ 



(8) 



fio ln(8/i/a) - l + k/2 ' 

Starting from a nearly potential configuration (small / 
value), if the current value could be progressively in- 
creased in function of time as in a classical electric cir- 
cuit, then a loss of equilibrium would occur when / would 



reach the maximum value of /cq. (/i). The circular cur- 
rent channel behaves in the same way as a straight cur- 
rent channel and its image current (FigurePc) a nd Sec- 
tion 13. 6p as proposed by Ivan Tend fc KuperusI ([1978). 
However, with an MHD evolution, the magnitude of the 
coronal current is rather determined by the magnetic 
field evolution, so its time evolution cannot be imposed a 
priori. Indeed, during a typical MHD evolution with im- 
posed photospheric velocities, the coronal cu rrent mag- 
nitude first grows then later decreases fe.g. lAlvl 119851 : 
iKlimchuk fc Sturrocklll989l: lAulanier et al.ll2005i r 

In previous studies, the magn etic field evolution is typ- 
ically assumed to be ideal (e.g. Ilsenberg fc Forbesll2007L 
and references therein) . With no magnetic flux emerging 
or canceling between the photospheric sources (located 
at X = iD), this implies the conservation of the mag- 
netic flux, F, through the area, S, defined between z = 
and the bottom of the current channel. F is given by 



II Sp^xdydz 



Lrl 



- A-Kcj) 1 - 



D 



^ih-a)^ + D^ 



(9) 



(10) 



with La = fj,Qh{\ii8h/a — 2) is the external inductance of 
the current channel (see Equation ([T])). The conservation 
of F, together with the time evolution of one parameter 
of the model (e.g. (j) or D) provides an evolution con- 
straint with the generic form /evoi.(^)- Its intersection 
with the equilibrium curve, /oq.(/i). determines the evo- 
lution of / in function of the evolving parameter (e.g. (j) 
or D). 

With a circular current channel, iLin et all ()2002l ) has 
shown that evolving D, with F preserved, does not lead 
to a loss of equilibrium, but rather to a self similar evolu- 
tion of the configuration (which is just a rcscalc in size). 
However, they also show that a decrease of <j> does lead 
to a loss of equilibrium. Figure EJa) is a graphical repre- 
sentation of this result. As in previous studies, the max- 
imum of LcqXh) defines the reference state with fluxes 
00 and _Fo. The evolution is parametrized by deflned 
by = /<;/) </'o I where is a decreasing function of time 
as o bserved after the emergin g phase of an active region 
(e.g. Ivan Driel-Gesztelvi et al . 2003) or before fllament 
eruption (e.g. ISchmieder et al.l 120081) . In order to have 
an equilibrium curve not evolving with (f> (to simplify 
the graphic), we draw /oq./0 (normalized to its maxi- 
mum value), as well as L^^oijcf). 

Starting the evolution with > 1, the constraint 
F = Fq has an intersection with the equilibrium curve 
before its maximum (e.g. at the point "a" for = 1.5 
in Figure |2ja)). This equilibrium is stable since the per- 
turbed equilibrium, with the constraint F = Fq, has a 
restoring force as shown with the arrows (see Section IITT] 
for the analysis of the perturbed equilibrium). 

As /0 decreases (due e.g. to magnetic flux cancelation 
with some photospheric flux brought from \x\ > D), the 
equilibrium height h increases. Contrary to what is ob- 
tained if the evolution of / is prescribed as in a circuit 
model, the conservation of magnetic flux in ideal MHD 
leads to the equilibrium to be still stable after the max- 
imum of Lf,q./4> (point "b"), up to the point "c". 

For larger h values, e.g. at point "a' " , the equilib- 
rimir is unstable to an ideal perturbation (i.e with F 
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Fig. 2. — The coronal electric current, /, divided by the photo- 
spheric magnetic flux 0, in function of the coronal height h, nor- 
malized to half the x distance, D, of the photospheric magnetic 
polarities. The equilibrium curve, Icq./(p{h), is in black, while the 
evolution constraint is shown in grey/dashed style for four values 
of the evolution factor /. Both / and Icq./(f> are normalized to 1 at 
the maximum of Icq./<j)- (^i^)) Case of a circular current channel 
(Figurc[TIb)), (c) case of a straight current channel (Figure [TJc) ) . 
(a,c) The photospheric flux, <j), is evolved with the factor f = f^. 
(b) The magnetic twist of the flux tube is evolved with the fac- 
tor / = /t • Arrows represent the direction of the force when the 
system is perturbed away from the equilibrium. 

preserved). However, this unstable region is not reach- 
able during the pre-cruptive evolution: the quasi-static 
evolution ends at point "c" with both an equilibrium be- 
coming unstable (see arrows in Figure [2l^a)) and without 
neighbor equilibrium, as is slightly more increased, so 
there is a loss of equilibrium. The evolution sequence 
with decreasing ends with a fast evolution driven by 
an outward force, so an eruption. It could be confined if 
a stable equilibrium would be present at greater height, 
e.g. due to the f ormation of a long curre nt sheet be- 
low the flux rope (|Forbes fc Isenberall991t) . or if a too 



strong magnetic tension of the overlying field would be 
present (as for the k ink-unstable eruption modeled by 
iTorok fc Klie^[2005h . 

3.3. Other Possible Evolutions on Intermediate Time 

Scales 

The analysis of the previous section supposes an ideal- 
MHD evolution to progressively shift the equilibrium to a 
point of the equilibrium curve where no neighbor equilib- 
rium exists (with the magnetic flux conservation). This 
is only one plausible scenario for solar eruptions (see Sec- 
tion [ij. Indeed, on the intermediate time-scale tb, mag- 
netic rcconncction can play a crucial role in transforming 
the coronal field. For example, this is the case with a 
progressive diffusion of photospheric magnetic polarities 
leading to reconnection at the PIL and the transforma- 
tion of sheared to twisted field lines, contributing to build 
up the twisted flux tube and the associated current chan- 
nel (as in the numeri cal simulations of Amari ct al. 2003; 
Fan fc GibsCT] | 2004t iMackav fc van BallegooiienI 12001 
Aulanier et all l2010t ). A progressive transformation of 



the coronal magnetic configuration is also expected due 
to reconnection of the current layers formed at separatri- 
ces a nd Quasi-Separatrix Layers (QSLs, e.g. iTitov et al.l 
I2OO2I ) during the intermediate time-scale tb- We con- 
clude that typically it is not obvious to justify an ideal- 
MHD evolution during the long pre-eruptive build-phase 
before an eruption, on time scales much longer than ta, 
and in particular using the evolution constraint of a pre- 
served flux below the flux rope (e.g. Equation PT!]) ). 

More generally, for a given observed configuration or 
in a 3D MHD simulation, it is difficult to determine the 
mai n evolutionary constra int, even with detailed analysis 
(e.g. I Aulanier et al.|[2010l ). This is not due to a lack of 
available data; rather, the difficulty is due to the com- 
plexities inherent in the 3D evolution of magnetic fields. 
This is illustrated by the formidable complexity of study- 
ing the loss of equilibrium even in a simplified configu- 

ration (e. g. a small bipole emerging in a bipolar field, 

iLin et al.l[200lh . 

3.4. Instability with a Circular Current Channel 

Let us illustrate another possible evolution than the 
ideal-MHD evolution on the time scale tb- This follow- 
ing case is selected mainly because of its simplicity to 
illustrate other possible evolutions. The current channel 
is associated with a twisted flux tube which has a finite 
twist due to the presence an axial magnetic field. The 
average coronal twist, T, is approximately related to the 
current / and the t oroidal flux, 0ti in the fl ux rope by 
(see Equation(9) in lTitov fc Demoulin|[l999l) : 

T = ti^Ih I (2,/.t) . (11) 

Let us suppose that the flux rope twist T is increasing, 
e.g. due to torsional Alfvcn waves coming from the con- 
vective zone, or as a consequence of reconnecting sheared 
loops (increasing the flux rope flux). Equation ([TT]) pro- 
vides a new evolution constraint, replacing the conserva- 
tion of F used in Section [321 As previously, the maxi- 
mum of loqXf^) defines the reference state with the flux 
(f>Q and the twist Tq. Here, we simply suppose that the 
photospheric held sources are not evolving so (j) = 4>o, but 
we still plot Icq. 1 4' for coherence with previous case. The 
evolution is parametrized by /t deflned by T = /t Tq . 
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Starting the evolution with /t < 1, the evolution con- 
straint, Equation ([TT|) . has an intersection with the equi- 
librium curve before its maximum (e.g. at the point "a" 
for /t = 0.5 in Figure [2i;b)). As in Section [33 this 
equilibrium is stable. As /t increases, the equilibrium 
height h increases, reaching point "b" then point "c" . If 
we only consider the evolution constraint, the evolution 
would reach point "d" , where no neighbor equilibrium is 
present when /t is further increased. As in the previous 
case (in Figure [^a)) a loss of equilibrium is present, but 
at larger height. 

However, the true physical evolution of the system is 
ending earlier, at point "c" , where the system is unstable 
on the short time scale ta with an ideal-MHD perturba- 
tion (preservation of the flux F in Equation (jlOp , while 
the evolution of T is negligible on the short time scale 
ta). 

3.5. Loss of Equilibrium or Instability? 

In the two previous evolutions (Figures [2ja) and (b)) 
the eruption occurs at the same location along the equi- 
librium curve, at the point "c" . However, they appear to 
have a different theoretical origin. In the first case (Fig- 
ure[2|a)) the progressive decrease of brings the system 
to point "c" where no neighbor equilibrium exits with a 
further decrease of f^. In the second case (Figure [2l^b)) 
the increase of /t brings also the system to point "c", 
but a neighbor equilibrium is present with a further in- 
crease of /t (up to point "d"). Simply, the equilibrium 
after point "c" is unstable with an ideal perturbation {F 
preserved), and as in the first case, an upward eruption 
is present after the system reaches point "c" . Indeed, on 
the time scale ta, the evolution of the two cases will be 
the same after they reach point "c" , since the same un- 
stable force (Equation ^) is acting with the constraint 
of an ideal-MHD evolution (so preserving F in Equa- 
tion (fTOj) ). Then, we argue that there is no point to dis- 
cuss whether there is a loss of equilibrium or a transition 
from a stable to an unstable equilibrium. 

More generally, for an observed coronal field evolu- 
tion, and even for an MHD simulation (where all physical 
quantities are available in the volume), we claim that, in 
most cases, it will be at least difficult, if not impossible, 
to define precisely the evolution constraint as in the pre- 
vious two cases. Indeed, with a slow enough driving, the 
system follows the equilibrium curve, and the extensions 
away from this curve, along the evolution constraint, are 
purely theoretical considerations which are available only 
if an analytical analysis is achievable. In a magnetic con- 
figuration which includes some separatrices or QSLs, so 
with some reconnection, we claim that an evolution con- 
straint cannot generically be constructed. However, a 
generic approach is to test the ideal-MHD stability all 
along the evolution (within the limit of numerical dissi- 
pation for MHD simulations). 

Because of the separation of the time scale of the coro- 
nal evolution (ta) from the longer times scales of the pho- 
tospheric evolution (tb, tc, see Section l3.ip . the mag- 
netic configuration cannot reach the equilibrium branch 
at an altitude h larger than that of the point "c" by 
a slow quasistatic evolution, which does correspond to 
observed pre-eruptive evolutions. Reaching a region of 
the equilibrium curve beyond "c" may be dynamically 
possible, but only for fast (e.g. Alfvenic) evolutions of 



magnetic configurations being out of equilibrium. This 
is, however, not observed in the Sun's atmosphere, even 
during flux emergence. Nevertheless, considering an an- 
alytical model of the equilibrium, one can always start 
a numerical simulation from any point along the equi- 
libriu m curve, stable or unstable (e.g. iTorok fc KliemI 
|2007|) . The unstable branch is unaccessible for a coro- 
nal field, except the vicinity of point "c" which can be 
reached with a small but finite velocity due to the slow 
evolution present in the pre-eruptive stage. Indeed, in- 
cluding such small velocity implies an evolution curve, 
h{t), which is in better agreement to observed h{t) in 
prominence eruptions, than letting the i nstability grows 
from an initial very small perturbation ([Schriiver et al.l 
[2008h . 

3.6. Loss of Equilibrium with a Straight Current 
Channel 

The evolution summarized Figures [^Ja) and (b) is ex- 
pected to be generic of magnetic configurations having 
at least one current channel which is not fully neutral- 
ized since the hoop for ce is generically prese nt with a 
curved current channel (iGarren fc Chenlll99l . Indeed, 
this has been shown in more complex configurations, 
even by including a complete photospheric line tying, 
i.e. not only fixing the "normal magnetogram" (with 
the inclusion of image current, see Section [2^ . but also 
fixing the photospheric p ositions of the current channel 
(jlsenberg fc Forbesl[2007l) . 

However, the above physical evolution is not limited 
to the presence of the hoop force. It is indeed generic 
of the Lorentz force created by any current channel. At 
a given position of the circuit, the magnetic force can 
be dominated cither by the magnetic field created by 
the current located in the vicinity of this position (hoop 
force) or by the c urrent at large distance (e.g. an image 
current ) . Indeed , Ivan Tend fc KuperusI (|1978D first pro- 
posed a catastrophe model of a straight current channel 
embedded in a potential field. 

As in Section [3?2l a potential field Bp is introduced to 
achieve an equilibrium. As previously, we select a bipolar 
field Bp created by two magnetic sources of fiux 4> located 
aX X = ±Z), z = 0, but now invariant by translation in 
the y direction. The x-component of Bp at the current 
location is: 

Bp^^^ -2(j)D{-K{h^ + D^))-^ (12) 

The equilibrium current, /oq. is given by / = in Equa- 
tion ®: 

Stt (j)Dh 

The equilibrium curve is closely similar to the one found 
for a circular current channel (compare panels a and c 
of Figure [2]). The main difference is that the maximum 
of Icq. is shifted to a larger height. This is mainly due 
to the different dependance with h of Bp ^ for a 2D and 
3D bipolc. Much closer equilibrium curves are obtained 
when the 2D bipole is replaced by a 2D quadrupole (giv- 
ing a potential- field dependance similar to Equation (|7|). 

The conservation of the magnetic flux below the cur- 
rent channel, per unit length along its axis (Equa- 
tion for the straight channel case is: 

^ UqI , 2h 26 ,h~a ,^ ^. 

F=—\n -tan"^ . (14) 

2tt a TT D 
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A major difference with the circular case is that the 
constraint of F conservation implies that / has a much 
weaker dependence on h (Figure[Uc)). It implies that the 
ideal-MHD instability occurs just after the maximum of 
^oq./0 (where = 1), when the photospheric polarities 
have weaken by only 4% (/^ = 0.96). Indeed, the ideal- 
MHD evolution leads to a non equilibrium at a location 
nearby to the one found by Ivan Tend fc KuperusI (|1978f ) 
with an evolution simply driven by an increase of the 
current /. 

As for the circular current channel, a loss of equilibrium 
at point "c" is only present if the magnetic evolution is 
fully ideal. More generally, whatever is the driver of the 
evolution on the time scale tb, the system is becoming 
unstable as it reaches point "c" , and it is ideally driven 
away from the equilibrium curve by the same force in the 
short time scale ta- 

3.7. Comparison to Previous Studies 

Based on our analysis above, we can now answer 
to the following question: is there a major differ- 
ence in the following approaches: "loss of equilib- 
rium" st udies with straight or circu l ar curr ent chan- 
nel (e.g. Ivan Tend fc Kuperuj Il978t iLin et a l. 20021, 
and the "torus instabilitv" (|Kliem fc Toro^OOd) pre- 
viou sly studied in t he tokamak laboratory experiment 
re.g. lBatema"a[T97ah ? 

First, the straight or circular current channels have 
very similar repulsion forces, r(/i)/^, implying the same 
kind of equilibrium curve. Indeed, for a circular current 
channel, rc{h) has a contribution from both the coro- 
nal part (z > 0) and from the image current (z < 0). 
Simply, both contributions can be combined in a single 
term. Equation (j4]), maskin g the contributiq i i of th e im- 
age current. More generallv lGarren fc ChenI (|1994[ ) have 
derived a general expression for r(h) for arbitrary current 
shapes, rc and Tg (Equations (|4|) and ([5])), are simply two 
limits of the same r(h), for circular and straight current 
channels^ 

Ncxt. lKliem fc Torokl ()2006l ) studied the "torus insta- 
bility" of a circular current channel imposing a constant 
current / or a constant flux F. The first case is direct ! 
comparable to the work of Ivan Tend fc Ku pcruj (jl97< 
and the second case to the work o f ILin et a l.'( 20021 ) The 
only significant difference is that iKliem fc Torokl ()2006l ) 
study the stability of the equilibrium curve (FigurejUJa)), 
but do not follow the evolution of the magnetic config- 
uration on intermediate time scale tb (so they cannot 
detect the presence or not of a loss of equilibrium). 

The main difference between the loss of equilibrium 
and the stability analysis is that, for the first case, a 
precise way how the system can evolve to instability is 
proposed, while, in the second case, one only tests the 
stability of a given equilibrium (which could be physi- 
cally inaccessible). However, with both type of analyses, 
an ideal instability is present at the same location of the 
equilibrium curve if the same equilibrium is analyzed. 
This last condition is not trivial if one allows the forma- 
tion of current sheets during the long-term evolution (on 
the time scale tb). 

4. EQUILIBRIUM STABILITY 

We study below the stability of the magnetic config- 
uration around an equilibrium position. The main as- 



sumption is that the magnetic force balance is the same 
at any position along the current channel, so that the 
magnetic force can written as in Equation (jS]). This in- 
cludes both circular and straight current channels in the 
same formalism. The stability of the configuration is first 
derived in this general framework in Section|431 with the 
constraint of magnetic flux conservation derived in Sec- 
tion 14.21 These results are used for the particular cases 
studied in previous section. Finally a parametric study 
of the stability is presented. 

4.1. Force Perturbation 

The magnetic force, f{a,h,I), on the current channel 
is described by Equation (j6|). A perturbation (da, d/i, dl) 
around the equilibrium, f{a,h,I) = 0, creates the force 
d/ as 



df 



da 



-da 



dh 



dh 



dl 
df 



d/. 



(15) 



During a perturbation the radius a is evolving as given 
by the internal force balance. With a linear force- free 
field inside the fiux rope and ideal MHD. lLin et al.l (|l998f ) 
found that a evolves as 1//. Here we include a more 
general variation, supposing a(/), and we introduce the 
index of variation of a with / 
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(16) 



so Tla 



1 for the internal evolution included in ILin et ahl 



(|1998[) . Then, Equation (|T5|) is rewritten as 
adf 



d/ 




I da 



df\ d/ 
'dl) dh 



dh. 



(17) 



pert. . 



where d//d/i|port. express how the current intensity is 
modified during the perturbation. 

If the perturbation is realized along the equilibrium 
curve, the force df has the same expression except that 
d//d/i|oq., computed along the equilibrium curve, re- 
places d//d/i|port in Equation pT)) . Also in this case, 
d/ = so 



dl ( _ adl 
dh^y'^^Ida 



a/\ d/ 
dl) dh 



0, 



(18) 



cq. 



where we suppose that the internal equilibrium has the 
same ria index. Using the equilibrium condition, / = 
rP ~ ~ 

rewritten as 



-Bp,x^ = 0, and Equation ([TH]), Equation ([TT]) is 



d/ = r / 1 - n. 



91ni 
^ d\x\( 







di 




) \ dh 


pert. 


~ dh 


eq./ 



dh. (19) 



The equilibrium is unstable when df and dh have the 
same sign. Then, an instability is present when the ab- 
solute value of the current in the perturbation decreases 
less rapidly with height than along the equilibrium curve. 
This is illustrated in Figure [2] 

On the short time-scale ta, the perturbation is de- 
scribed by ideal MHD, with a preservation of the flux 
distribution at the boundary z = (so 4> =constant in 
the examples of Section 13.21 and 13. 6p . Then, a pertur- 
bation corresponds to a small excursion away from the 
equilibrium curve with the constraint of magnetic flux 
conservation in Figures [2ja) and (c) . The tangent point 
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"c" between the equilibrium curve and the constraint of 
flux conservation define both the limit of the stable re- 
gion and the non equilibrium point (in an ideal evolution 
on the time-scale tb). 

Using the equilibrium condition to specify d//d/i|cq., 
the force perturbation, Equation (|17p . is rewritten as 



dh 



rl 1 



d In r 
d\na 



dl 
dh 



pert. 



dh 



We define the decay index of the potential field as 

dln\B„ 

nBp = - 



, dr 
dh ■ 
(20) 



d\nh 



(21) 



This decay index is introduced in Equation (PH)) by di- 
viding df by Bp ^^I, or equivalently by —rl^, and by 
multiplying by h. Finally, the equilibrium is unstable 
(i.e. df/dh > 0) if 

"Bp > ?^Bp,crit. ="-r + m 

dlnr 



dlnh 



1 



51nr\ dlnl/l 
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(22) 

pert. 



SO if the potential field decreases fast enough with height. 
The instability threshold depends on how fast the repul- 
sion decreases with height (given by the decay index n,-) 
and on how much the current is allowed to decrease dur- 
ing the perturbation (given by the decay index ni). 

4.2. Constraint of Flux Conservation 

The decay index ni is computed from an ideal MHD 
constraint, as follows. On the short time-scale r^, the 
perturbation is described by ideal MHD, with the preser- 
vation of the coronal magnetic flux F present below the 
flux rope. With a small perturbation dh, Equation ^ 
implies 



dF : 



die 
dh 



L. 



d/ 
dh 



P 



, (23) 



pert. 



where P is the perimeter (or length) of the full current 
channel, including its image. 

The external inductance (computed with = 0) de- 
pends generically on the spatial extension of the current 
channel (described here only by h) and on the thickness 
of the current channel (its radius a). Then, the variation 
of the external inductance with height is 
dLc dLc dLc da dl 
^~dh^ 



dh 



dh 



da dl dh 
dlnL, 



pert. 

din I 
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(24) 



pert. 



The conservation of the magnetic flux, Equation ((23)) . 
together with Equation ([24]) imply the following decay 
index 



1 



ni = 



In r 
^ d\na 



1 — 71; 
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dlnh 



ie 



(25) 



The last factor in the previous equation is due to the con- 
servation of the magnetic flux from the current / minus 
the flux from the potential field passing below the current 
channel. The fraction includes the effect of the channel 
expansion on the repulsion function at the numerator and 
on the conservation of flux at the denominator. They are 
both larger than 1 for ria > 0, and they arc a growing 
function of Tia- Then, the dependance of r and Lg on a 
have opposite effects on the stability. 



4.3. Example of a Circular Current Channel 

We apply in this subsection the above results to a cir- 
cular current channel (Figure [ijb)). The external induc- 
tance La is given by Equation ([1} , without the contribu- 
tion of the internal inductance (/i ~ 0). The repulsion 
r is given by Equation (|4]) and P = 27r/i. The contribu- 
tion from the potential fleld to the decay index nj is (for 
ria = 0, i.e. for a constant during the evolution): 

_ 1 / l + k/2 \ 

" 2 \^ ln{8h/a)~2j ' 

The contribution from the fleld created by the current 
channel to the decay index ni is (for = 0) 



dhiLc 
dlnh 



1 



" ^ + ln{8h/a) - 2 ' ^^^^ 

so exactly twice the negative contribution from the po- 
tential fleld for If = 0. For a h, i.e. a very thin 
current channel, the denominator, ln(8/i/a) — 2, is rela- 
tively large so that the above contributions are close to 
— 1/2 and 4-1, respectively. Including Equations (|26l27p 
in Equation and using Equation ([1]), the decay in- 
dex Hi for the circular current channel is 

^ iu + k/2 + n,){u-k/2) 

2iu~l + n,){u + lf/2) ' 

with the notation 

u = ln{8h/a) - 1 . (29) 
Using Equation ([5]), the decay index of the repulsion r 

dlnr 



is: 



dlnh 



1 



1 



(30) 



u + k/2 

From Equation (P^ . is the critical decay index of Bp 
for instability if the current / would be preserved during 
the perturbation. Ur is always lower than 1, especially 
for flux rope with large radius a. 

Combining the results of Equations p8|30p , the insta- 
bility condition for a circular current channel is 



"Bp > 2 



{l + k/2){u-2 + 2n^ + k/2) 



(31) 



2{u + k/2){u-l + n^) 

In the limit of a very thin current channel (a <C h), the 
instability th reshold is close to 3/2, as found in tokamak 
studies (e.g. iBatemanI [19781 ). This cor responds to the 
"toru s instability" for solar eruptions (jKliem fc Torcikl 
|2006[ ). but with a different correction term to 3/2, as we 
have not supposed a self similar expansion of the current 
channel, but rather a dependance a(/). 

4.4. Example of a Straight Current Channel 

We follow the same derivation than in previous sub- 
section but for a current channel formed by two parallel 
lines (Figure [ijc)). In this geometry, the external induc- 
tance Le is given by Equation again with Zj = 0. The 
repulsion r is given by Equation ([5]) and P = 2Ay (with 
Ay ^ h,D). For ria = 0, the contribution from the fleld 
created by the current channel to the decay index ni, 
Equation (j25[) , is exactly twice the contribution from the 
potential field, but with opposite sign (as above for the 
circular channel, with l{ = 0) 

^ (32) 
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Fig. 3. — Dependence of the decay indexes nr 
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n\ = — 9 in In /i and nep = — 9 In [Bp^xl/fln /i for a straight 
(resp. circular) current channel in the left (rcsp. right) column. 
The critical index riBp.crit. = n-c + n\ for instability is shown 
with a black continuous line. The decay indexes of 2D and 3D 
potential bipolcs are shown in both columns for reference. The 
dotted vertical line indicates the critical height (point "c" in Fig- 
ure [2]l. All panels are drawn with the normalized equilibrium ra- 
dius a/D = 0.1 and for h > a. The top (resp. bottom) panels have 
Ha = —dlna/dln \I\ = (resp. na = 1). 
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Fig. 4. — The critical index riBp.crit. = rii + ni for instability 
with various equilibrium radius a normalized by D (half distance 
between photosphoric field concentrations). The results are for a 
straight (resp. circular) current channel in the left (resp. right) col- 
umn. The flux rope expands more during the perturbation (larger 
ria value) from top (ria = 0) to bottom (ria = 1). The nBp,crit. 
curves are drawn for h > a. The decay indexes, riBp of 2D and 3D 
potential bipoles are shown in both columns for reference. 

With Equation ([Sj, the decay index of the repulsion r 
is simply = 1. Combining the above results, the in- 
stability condition. Equation ([22]). for a straight current 



"Bp > 1 



1 



2(ln(2/i/a) + na 



(33) 



In the limit of a very thin current channel (a ^ /i), 
the instability thres h old i s close to 1, as found by 
Ivan Tend fc KuperusI (|1978l ). 

4.5. Comparaison of Circular and Straight Current 
Channels 



From Equation ([22)) . the critical decay index of the po- 
tential field for instability, nep.crit. has two contributions: 
the decay indices of the repulsion, n,-, and the decay in- 
dex of the current during the ideal perturbation, nj. The 
main difference of stability between the above two cur- 
rent channels is a much lower index ni for a straight 
channel (Figure ^ . This is due to the low dependance 
of Lc on the height h (compare Equation ^ to Equa- 
tion (H])). In the limit a <^ h, this is the origin of a more 
stable circular current channel (nBp,crit. ~ 3/2) than a 
straight current channel (?T.Bp,crit. ~ !)• However, this 
limit is not applicable to the eruptive coronal configura- 
tions since it requires extremely thin current channels: 
for example, even with a/D ~ 10"'^, this limit is only 
weakly approached (Figure^]). This is due to the ln(/i/a) 
dependance in both Equations ((3T|) and (|33p . 

Moreover, the above difference is partly compensated 
by a lower decay index of the repulsion, rii-, for a circular 
channel (as the repulsion, r(/i), is decreasing slower with 
height, especially for large values of a). The net result 
is that the critical decay index, nBp,crit., for circular and 
straight channels, with finite radius a, are much closer 
than in the limit a h. This is illustrated in Figure [3] 
with a/D = 0.1, which still corresponds to a relatively 
narrow channel. This effect is amplified with a flux rope 
having a larger expansion during the perturbation, so a 
larger ria value, because the contribution of ni to nBp,crit. 
is reduced. Indeed, even with a relatively thin flux rope, 
a/D = 0. 1 , and an internal expansion rate ria = 1, as in 
iLin et all (|l998f ). JT-Bp.crit. values are close for circular and 
straight current channels: w 1.3 and ~ 1.1, respectively 
(Figure [31Jc) and (d)). They correspond to a comparable 
critical height {h/D « 1) for straight and circular chan- 
nels with a 2D and 3D bipole (Equations ([7]) and ([T^ . 
respectively) . 

In previous studies, the instability threshold 7iBp,crit. 
was typical ly taken in the limit of very t h in current chan- 
nels (e.g. Ivan Tend fc KuperusI 119781: iKliem k T6m^ 
|2006[ ). The above analytical theory is indeed done in 
the limit of thin current channels (typically a/D < 
0.1). In fact, relatively broad current channels are ex- 
pected in the coronal with magnetic extrapolations (e.g . 
ISchriiver et al .1120081: [Savcheva fc van Ballegooiienll2009l) 
and are present in MHD simulations (e.g. iFan &: GibsonI 
I200llm5k fc Klienill2007l:IAulanier et al.ll2010D . So. we 
also show the approximative results for a broad channel 
case, a/D = 0.5. For a straight current channel, the flux 
conservation provides an increasing stabilizing effect as 
the channel radius, a, is increasing (Figure^]). Indeed, for 
a broad straight channel {a/D > 0.1), ?T.Bp,crit. can reach 
a value comparable, or even larger in some casesthan the 
one obtained for the corresponding circular channel with 
the same parameter values (Figure jH^e) and (f)). We 
conclude that the circular and straight current channels 
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have typically comparable instability threshold values for 
the range of parameters expected in the corona. 

The internal evolution of the current channel during 
the perturbation (i.e. the effect of ria), has a common ef- 
fect on the stability of a circular and of a straight current 
channel because the stabilizing term provided by the flux 
conservation is generally decreased with an increase of ria 
(Figures I3I4|) . However, changing the channel thickness 
(a/D) has an opposite effect on the instability threshold 
for a circular and for a straight current channel (compare 
the panels in Figure |4] where the a/D curved are ordered 
oppositely in the two columns). For a straight channel, 
rir = 1, so riBp^crit. = + Til is affected only by the de- 
pendance of ni on a/D. But, for a circular channel, the 
n,- depcndancc on a/D is important, and dominates the 
contribution of ni, as soon as ria is slightly positive. It 
implies a dependance of riBp.crit. on a/D for a circular 
channel that is opposite to that for a straight channel. 

5. CONCLUSION 

How to destabilize a coronal magnetic configuration is 
a key-issue of CME research. Among several possibilities, 
two candidates are a loss of equilibrium and a torus insta- 
bility occuring during the evolution of the magnetic con- 
figuration. Both have been initially developed with the 
approximation that the coronal currents are restricted to 
a non-neutralized current channel, and both theory were 
further analyzed with MHD simulations, relaxing part 
of the initial approximations of the analytical develop- 
ments, but at the expense of not covering the parameter 
space. 

In this study, we revisit both analytical theories and 
compare their approaches for the two simple configura- 
tions where their results apparently differ: a straight 
and a circular current channel. A loss of equilibrium 
is typically, but not always, present in both configura- 
tions when an ideal-MHD evolution is imposed during 
the long-term evolution of the magnetic configuration. 
However, when a loss of equilibrium occurs, the magnetic 
configuration is also ideally unstable. From the results of 
Sections|3]and|4l wc conclude that both approaches are in 
fact both compatible and complementary. In particular 
they agree on the position of instability, if no significant 
current sheets are formed during the long-term evolution 
of the magnetic configuration. Moreover, slow resistive 
processes, e.g. tether cutting, are probably occurring all 
the way before an eruption occurs. Therefore, wc con- 
clude that the analytical theory is most useful in deriving 
an instability threshold with the constraint of ideal MHD 
evolution on a short time scale (coronal Alfven time) for 
a given magnetic equilibrium. 

Wc also compare the physical origin of the instability of 
straight and circular current channels. In order to model 
the negligible evolution of the magnetic flux crossing the 
photosphere on the coronal Alfven time scale, a theoret- 
ical image current is introduced below the photosphere 
(Figure [T]). For a straight current channel, the repul- 
sion of the image is balancing the Laplace force between 
the coronal current and the potential field (associated 
to the photospheric field distribution). For a circular 
current channel, the repulsion of the nearby coronal cur- 
rent is also present (called hoop force). However, since 
the repulsion force depends only on the global curva- 
ture radius and on the thickness of the current chan- 



nel for a circular channel, it could lead the false conclu- 
sion that the repulsion force has a different origin for the 
str aight and circular curr ent channels. In fact, as shown 
bv lGarren fc Chenl ()1994D . both the coronal and the im- 
age currents generically contribute to the repulsion force 
of a current channel. Both terms actually combine in a 
single expression for a circular channel, while there is no 
contribution of the coronal current for a strictly straight 
current channel. The circular and straight current chan- 
nels are simply two limits of the general case with specific 
properties. 

The instability occurs when the potential magnetic 
field decreases fast enough with height, more precisely 
when its decay index, riBp as defined by Equation (|2ip . 
is larger than a critical value riBp.crit.- At the limit of 
extremely thin current channels nBp,crit. = 1 and 1.5 for 
a straight and circular current channel, respectively. In 
fact, we show that this difference is not due to the differ- 
ence in the repulsion force, but that it is due to the con- 
straint of ideal MHD (conservation of the coronal mag- 
netic flux below the current channel). Moreover, with 
the sole contribution of repulsion force to the instability 
threshold (i.e. ni = in Equation ([22| ). first nBp,crit. < 1 
for a circular current channel while nBp,crit. = 1 for a 
straight channel, and second riBp.crit. approaches 1 for 
both circular and straight channel as the channel be- 
comes very thin. This further indicates that there is no 
real difference in the origin of the repulsion force for a 
straight and circular current channels. 

We conclude that the same physics is involved in 
the instability of circular and straight current channels, 
and that they are just two particular limiting cases of 
more general current paths. For the typical range of 
current-channel thickness expected in the coronal, and 
present in MHD simulations, and for a current chan- 
nel expanding during an upward perturbation, nBp,crit. 
has close values for both circular and straight current 
channels (in the range [1.1,1.3], Figure Ul^e) and (f)). If 
the current channel would not expand, the decay in- 
dex riBp.crit. would be higher, typically in the range 
[1.2,1.5], but still not so different in both cases (Fig- 
ure m^a) and (b)). Similar critical indexes have been 
found in MHD simulations starting from a initial equi- 
libriur n having a coronal current channel close to half 
torus (iTorok fc Klieml[2007t ISchriiver et al.l[2008l) . Oth- 
erwise, from the measurement of the height of a set of 
quiescent pro minences, combuied wi th potential field ex- 
trapolations, |FiiiEE2^SSi3 (|2001[ ) found nBp,crit. ~ 1- 
This threshold is closer to riBp^crit. of the straight current 
channel as expected since quiescent prominences are hor- 
izontally extended structures. 

In an MHD simulation, where a fiux rope and its asso- 
ciate current channel is progressively formed due to pho- 
tospheri c motions, flux cancel ation, and magnetic recon- 
nection, lAulanier et al.l (|2010( ) found an unstable config- 
uration when TT-Bp.crit. ~ 1-5. There, the fiux rope height 
satisfied h/D w 1. This was in favor of the "torus insta- 
bility". However, with the above results, this threshold 
would require that the current channel is almost rigid 
during the perturbation (i.e. na = as in Figure IUb)). 
It is not obvious that this condition is met in a low-/? 
magnetic field. More relevant is probably the role of the 
anchorage of the current channel at fixed photospheric 
positions during the stability analysis, a constraint not 
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present for both the straight and circular models stud- 
ied above, but included in the non-equilibrium study of 
llsenberg fc Forbei (j2007() . The present analytical theory 
is also over simplifying the coronal current distribution to 
only one current channel, while at least partial current 
neutralization as well as other current layers arc typi- 
cally present in MHD simulations. This may also raise 
the critical index, riBp.crit., to larger values (e.g. nBp,crit. 
as high as 1.9 was found in the MH D simulations of flux 
emergence bv lFan fc GibsM] [20071 ) . The precise under- 
standing of the instability threshold is important for de- 
termining when a CME would occur. This will be the 
object of further developments of the analytical theory. 

It would also be desirable to derive the critical index 
from observations of eruptive prominences and sigmoids, 
the pre-cruptive altitudes of which can either be mea- 
sured using two vantage points (e.g. using the pair of 
STEREO imagers) or when they cross the solar limb. 



At first approximation, the background coronal magnetic 
field would then have to be extrapolated in the potential 
field approximation, using photospheric magnetograms 
either taken on the same day of the eruption if possible, 
or a few days before if no magnetogram is available. Such 
a survey of various erup tive solar features wou ld extend 
the work carried out by iFilippov fc DenI ()200lD , who fo- 
cused on long and high altitude quiescent prominences, 
that mostly concern the straight channel model. 
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